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We investigate the vacuum expectation value of the energy-momentum tensor associated 
with a massive fermionic field obeying the MIT bag boundary condition on a spherical 
shell in the global monopole spacetime. The asymptotic behavior of the vacuum densities 
is investigated near the sphere center and surface, and at large distances from the sphere. 
In the limit of strong gravitational field corresponding to small values of the parame- 
ter describing the solid angle deficit in global monopole geometry, the sphere-induced 
expectation values are exponentially suppressed. 
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1. Introduction 

Global monopoles are spherically symmetric topological defects created due to 
phase transition when a global symmetry is spontaneously broken and they have 
important role in the cosmology and astrophysics. The simplest theoretical model 
which provides global monopoles has been proposed by Barriola and Vilenkin. ^ 
Neglecting the small size of the monopole's core, the corresponding metric tensor 
can be approximately given by the line element 



where the parameter a^, smaller than unity, depends on the symmetry breaking 
energy scale and codifies the presence of the global monopole. This spacetime cor- 
responds to an idealized point-like global monopole. Here we shall calculate the 
vacuum expectation values of the energy-momentum tensor for fermionic fields 
obeying MIT bag boundary condition on the spherical shell concentric with the 
global monopole. 
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2. Vacuum Energy-Momentum Tensor Inside a Spherical Shell 

The dynamics of a massive spinor field on a curved spacetime is described by the 
Dirac equation 

»7'^(5^+r^)^-MV' = , (2) 

where 7^ are the Dirac matrices and — \^u^ fil'^ is the spin connection with 
being the standard covariant derivative operator. Our interest in this paper 
will be the vacuum expectation values of the energy-momentum tensor induced by 
a spherical shell in the global monopole spacetime. We shall assume that on the 
sphere surface the field satisfies bag boundary condition: 

(1 + i7^nb,0 V = , r = a, (3) 

where a is the sphere radius, n;,^ is the outward-pointing normal to the boundary. 
Expanding the field operator in terms of the complete set of single-particle states 
ijj^p ^1 and making use the standard anticommutation relations, for the vac- 
uum expectation values of the energy-momentum tensor one finds the following 
mode-sum formula 

(0|T^.|0) = 5]r^.{4-)(x),v^(-)(x)} , (4) 

where |0) is the amplitude for the corresponding vacuum. For the geometry un- 
der consideration the eigenfunctions are specified by the set of quantum numbers 
(3 = {kjma), where j — 1/2, 3/2, . . . determines the value of the total angular mo- 
mentum, m = —j, ■ ■ ■ ,j is its projection, and cr = 0, 1 corresponds to two types of 
eigenfunctions with different parities. These functions have the form 



u; = ±E, E=Vk^ + M^ , = 1±1/? - !^ (6) 

a 2 

where n = f/r, a ~ (cr^, cr^, ct'^) with the curved space Pauli 2x2 matrices . 
In Eq. |SJ| Zjj(x) represents a cylindrical function of the order v and rij;,„ are the 
standard spinor spherical harmonics^ with I — j ~ n„/2. 

For the region inside a spherical shell one has Z^, (x) — J^, (x) , where Ji, (x) is the 
Bessel function. The imposition of the boundary condition on the eigenfunctions 
||SJ) leads to the following equations for the eigenvalues 

ljka)^0, (7) 

where for a given function F(z) we use the notation 

F{z) EE zF'{z) + (a^ + s^Vz^ + A^' - (-1)"'^) F{z) , a = 0, 1, (8) 

with Suj = sgn{uj) and /i ~ Ma. Let us denote by Xi^^.s = ka, s — 1,2, ... , the 
roots to equation Q in the right half plane, arranged in ascending order. By using 
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the standard integral for the Bessel functions, for the normahzation coefficient one 
finds 



z Vz^ + a2Af2 + aM 
2a2a2 + a2Af 2 



with 



(9) 



(10) 



Because the spacetime is spherically symmetric and static it follows that the 
vacuum energy-momentum tensor is diagonal: 



(0 |r;| 0) = diag(£, -p, -Pi, -p^) , 



(11) 



with the energy density e, radial, p, and azimuthal, p_L, pressures being func- 
tions on the radial coordinate only. As a consequence of the continuity equation 
V 1/ (O |r^| O) = 0, these functions are related by the equation 



r^+2(p-pi) = , 



(12) 



which means that the radial dependence of the radial pressure necessarily leads to 
the anisotropy in the vacuum stresses. Below we will give formulas for e and p. The 
azimuthal pressure is expressed via the radial one by using formula H12|l . 

Substituting eigenfunctions © into Eq. the summation over the quantum 
number m can be done by using standard summation formula for the spherical 
harmonics. For the energy-momentum tensor components one finds 



q{r) 



4 



i=i/2 



o-=0,l s=l 

g = £, 

where we have introduced the notations 

/(t' [z. My)\ = z \{^/z^T^ - fi)jUy) + iVz^T^ + lAJl+n^ (y) 



f!^'J[z,My)] 



^z^ + Ai2 



y . 



(13) 

(14) 
.(15) 



Applying to the sums over s the generalized Abel-Plana summation formulJSI, 
the components of the vacuum energy-momentum tensor can be presented in the 
form 



qir) = Qmir) + qb{r), 



(16) 



where Qmif) does not depend on the radius of the sphere a and is the contribution 
due to unbounded global monopole spacetime. The corresponding quantities for the 
massless case are investigated in Ref. ^ The second term on the right of formula 
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lfTB|) is induced by the presence of the spherical shell and can be presented in the 
form 



(17) 



with 



FiP) [x,Uy)] = - - . (19) 

y 

Here and below I ~ j + 1/2, 1^ = 1/1= //a + 1/2, and for given functions f{x) and 
g{x) we use the notation 

W [f{x),g{x)] = [xf{x) + + iy)f{x)] [xg'{x) + (/i + iy)g{x)] 

+ [x^-^^)f[x)g{x). (20) 

It can be easily checked that for a massless spinor field the boundary-induced 
part of the vacuum energy-momentum tensor is traceless and the trace anomalies 
are contained only in the purely global monopole part without boundaries. 

Now we turn to the consideration of various limiting cases of the expressions 
for the sphere-induced vacuum expectation values. In the limit r 0, for the 
boundary parts p7|l the summands with a given I behave as r^'/"^^, and the 
leading contributions come from the lowest I — 1 terms. Making use standard 
formulae for the Bessel modified functions for small values of the argument, for the 
sphere-induced parts near the center, r <C a, one finds 

^ .-^«-^(r/2«)^-^ r °^ ^ ^^^ W[L{x\KAx)]-^, 



TT 



-2a-4(r/2a)°"" xi+^dx W [h{x), K,{x)] 



^ 2a(2 + a)r2 (i i) ,/x^ W [hixihi^)] ' ^^^^ 

where v ~ l/a + 1/2 and V{x) is the gamma function. Hence, at the sphere center 
the boundary parts vanish for the global monopole spacetime (a < 1) and are finite 
for the Minkowski spacetime {a = 1). Note that in the large mass limit, /i ^ 1, the 
integrals in Eqs. f l21() . H22() are exponentially suppressed by the factor e^^^. 

The boundary induced parts of the vacuum energy-momentum tensor compo- 
nents diverge at the sphere surface, r — > a. In order to find the leading terms of the 
corresponding asymptotic expansion in powers of the distance from the sphere sur- 
face, we note that in the limit r ^ a the sum over I in H17|l diverges and, hence, for 
small 1 — r/a the main contribution comes from the large values of I. Consequently, 
rescaling the integration variable x vx and making use the uniform asymptotic 
expansions for the modified Bessel functions for large values of the order, to the 
leading order one finds 

_ /x+1/5 _ 1/5 - 2m 

\2'K'^aya — ry 247r^a^(a — r)^ 
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Notice that the terms in these expansions diverging as the inverse fourth power of 
the distance have canceled out. This is a consequence of the conformal invariance 
of the massless fermionic field. It is of interest to note that the leading terms do 
not depend on the parameter a and, hence, are the same for the Minkowski and 
global monopole bulks. In Fig. ^ we have presented the dependence of the Casimir 
densities, a'^qb{r) on the rescaled radial coordinate r/a for a massless spinor field 
on the Minkowski bulk (left panel) and on the global monopole background with 
the solid angle deficit parameter a = 0.5 (right panel). The vacuum energy density 
and pressures are negative inside the sphere. 

0.1 ^ ^ ^ ^ ^ 0.1 I ^ ^ ^ 1 ry • • 1 




Fig. 1. Vacuum energy density, a*e (curve a), azimuthal pressure a*px (curve b), and radial 
pressure a*p (curve c) for a massless spinor as functions on the ratio r/a inside and outside 
a spherical shell in the Minkowski spacetime (a = 1, left panel) and for the global monopole 
spacetime with a = 0.5 (right panel). 



Now let US consider the limit a ^ 1 for a fixed value r < a. This limit cor- 
responds to strong gravitational fields. In this case one has v f^i l/a ^ 1. The 
main contribution to the sum over I comes from the summands with I — 1 and 
the boundary parts of the vacuum energy-momentum tensor components behave as 
exp[— 2 ln(a/r)/a] with pb/p±b ~ oc. Hence, for a ^ 1 the boundary-induced vac- 
uum expectation values are exponentially suppressed and the corresponding vacuum 
stresses are strongly anisotropic. Fig.|21shows that the nonzero mass can essentially 
change the behavior of the vacuum energy-momentum tensor components. In this 
figure we have depicted the dependence of the boundary induced quantities a'^qb{r) 
on the parameter Ma for the radial coordinate r = 0.5a. The left panel corresponds 
to the sphere in the Minkowski spacetime {a = 1) and for the right panel a — 0.5. 

3. Vacuum Expectation Values Outside a Spherical Shell 

Now let us consider the expectation values of the energy-momentum tensor in the 
region outside a spherical shell, r > a. The corresponding eigenfunctions have the 
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Fig. 2. Boundary-induced vacuum expectation values a'^qi,{r), q = e,p,pj_, as functions on = 
Ma for r/a = 0.5. The curves a, b, c correspond to the energy density (e), azimuthal pressure 
{p±), and radial pressure (p ), respectively. For the left panel o = 1 (Minkowski spacetime) and 
for the right panel a = 0.5. 



form ^ with 

Z,{kr) 



gp{ka, kr) = J^{kr)Y,j{ka) — Y^{kr)J^{ka) , 



(24) 



where Y^{z) is the Neumann function, and the functions with tilda are defined 
as (IHl. Now the spectrum for the quantum number k is continuous and for the 
normaUzation coefficient we obtain 

^2 ^ k{u: + M) 

" 2a'^Lu[J^^{ka)+Yl{ka)] ' ^ ' 

Substituting the eigenfunctions Q into the mode-sum formula and taking into 
account Eqs. (|24|l and H25() . we can see that the vacuum energy-momentum tensor 
has the form The diagonal components are determined by formulae 



q{r) 



i=i/2 



(2j 



^=0,1-^0 



' ^ fall [x,gu^ix,xr/a)] 
JUx)+Ylix) ' 



p, (26) 



where the expressions for faf!^ [x,g^^{x,xr/a)] are obtained from formulae H14|l . 
((T31) by replacements 



(27) 



To find the parts in the vacuum expectation values of the energy-momentum tensor 
induced by the presence of the sphere we subtract the corresponding components 
for the monopole bulk without boundaries. After rotation of the integration contour 
and introducing the Besscl modified functions, for the boundary-induced parts one 
obtains 

1 v"^ , /"°° x^dx W \I^(x), Kiy{x)] ,„) , , , 
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Here the expressions for the functions Fal) [x,K^(y)] are obtained from formulae 
(fTH|l . ((T^ by replacements Iu{y) K^{y) and I,^-i{y) —Ky-i{y). As for the 
interior region, it can be seen that in the limit of strong gravitational fields, a <C 1, 
the boundary induced vacuum expectation values are exponentially suppressed by 
the factor exp[— (2/q;) ln(r/a)], an the corresponding vacuum stresses are strongly 
anisotropic: Pb/pi_b ^ a. 

In Fig. n we have plotted the dependence of the vacuum energy density and 
stresses on the radial coordinate for a massless spinor field outside a sphere on the 
Minkowski bulk (left panel) and on background of the global monopole spacetime 
with a = 0.5. As seen form these figures, the energy density and azimuthal pressure 
are positive outside a sphere, and the radial pressure is negative. The latter has the 
same sign as for the interior region. 

For the case of a massless spinor the asymptotic behavior of boundary part H28|l 
at large distances from the sphere can be obtained by introducing a new integration 
variable y = xr /a and expanding the subintegrands in terms of a/r. The leading 
contribution for the summands with a given I has an order (a/r)^""'"^ and the main 
contribution comes from the I = 1 term. The leading terms for the asymptotic 
expansions over a/r can be presented in the form 

where 

^^ = 4^;^' U = -T^^ /p^ = i- (30) 

As for the interior components, the quantities 1)28(1 diverge at the sphere surface 
r = a. Near the surface the dominant contributions come from modes with large 
I and by making use the uniform asymptotic expansions for the Bessel modified 
functions, the asymptotic expansions can be derived in powers of the distance from 
the sphere. The leading terms of these asymptotic expansions are determined by 
formulae 

/ \ 1/5 — 5/i , , 1/5 — 2/i 

^bi^) ^ ir, 2 / ^ ' Pb^^) OA 2 2i YI ' ^1) 

Recall that near the sphere the interior energy density is always negative. As we 
see, the leading terms for the radial pressure are the same for the regions outside 
and inside the sphere. For the azimuthal pressure these terms have opposite signs. 
In the case of the massless spinor field the same is true for the energy density. 



4. Concluding Remarks 

In this paper we have analyzed the fermionic Casimir densities induced by a spheri- 
cal shell in a idealized point-like global monopole spacetime. The boundary-induced 
expectation values for the components of the energy- momentum tensor are given by 
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formulae (1171 ) and (|28|l for interior and exterior regions, respectively. These expres- 
sions diverge in a non-integrable manner as the boundary is approached. The energy 
density and azimuthal pressure vary, to leading order, as the inverse cube of the 
distance from the sphere, and near the sphere the azimuthal pressure has opposite 
signs for the interior and exterior regions. For a massless spinor the same is true for 
the energy density. The radial pressure varies as the inverse square of the distance 
and near the sphere has the same sign for exterior and interior regions. The leading 
terms of the corresponding asymptotic expansions near the sphere do not depend 
on the solid angle deficit parameter and are the same for these two cases. Near the 
sphere the interior energy density is negative for all values of the mass, while the 
exterior energy density is positive for Ma < 0.04 and is negative for Ma > 0.04. 
Near the sphere center the dominant contributions come from modes with I — 1 
and the sphere-induced vacuum expectation values vanish for the global monopole 
spacetime and are finite for the Minkowski bulk. At large distances from the sphere 
the components of the vacuum energy-momentum tensor go to zero as (a/r)^/"^^. 
In the limit of strong gravitational field, corresponding to small values of the param- 
eter a, describing the solid angle deficit, the boundary-induced part of the vacuum 
energy-momentum tensor is strongly suppressed by the factor exp[— (2/q;)| ln(r/a)|] 
and the corresponding vacuum stresses are strongly anisotropic: pi, ^ ap±h. Note 
that this suppression effect also takes place in the scalar case. 
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